#open(BUFF,"fwd_ambi");
open(BUFF,"rev_ambi");
open(IMP,"impu_ambi");
while(<BUFF>)
{
	chomp($_);
	$imp=<IMP>;
	chomp($imp);
	@a=split(" ",$_);
	@imp=split(" ",$imp);
	$snp=shift(@a);
	shift(@imp);
	$num=0;
	for($i=0;$i<@a;$i++)
	{
		$a1=$a[$i];
		$b1=$imp[$i++];	
		$a2=$a[$i];
		$b2=$imp[$i++];
		$a3=$a[$i];
		$b3=$imp[$i];
		if($a1 >= 0.5)
		{
			$val1=1;
		}
		elsif($a2 >= 0.5)
		{
			$val1=2;
		}
		elsif($a3 >= 0.5)
		{
			$val1=3;
		}
		else
		{
			$val1=4;
		}
		if($b1 >= 0.5)
                {
                        $val2=1;
                }
                elsif($b2 >= 0.5)
                {
                        $val2=2;
                }
                elsif($b3 >= 0.5)
                {
                        $val2=3;
                }
                else
                {
                        $val2=4;
                }
		if($val1 != $val2 )
		{
			$num++;	
		}
		#$num=$num/513;
		#print "$snp $num\n";	
		#print "$val1 $val2\n";
	}
	$num=$num/513;
	print "$snp $num\n";
}
